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Abstract. Discretizations of the Bogoyavlensky lattices are introduced, be- 
longing to the same hierarchies as the continuous-time systems. The construction 
exemplifies the general scheme for integrable discretization of systems on Lie alge- 
bras with r-matrix Poisson brackets. An initial value problem for the difference 
equations is solved in terms of a factorization problem in a group. Interpolating 
Hamiltonian flow is found. 



1 Introduction 



The subject of integrable symplectic maps received in the recent years a 
considerable attention. Given an integrable system of ordinary differential 
equations with such attributes as Lax pair, r-matrix and so on, one would like 
to construct its difference approximation, desirably also with a (discrete-time 
analog of) Lax pair, r-matrix etc. Recent years brought us several successful 
examples of such a construction [1-10]. 

Recently, stimulated by the results of [^], there was formulated a 
general recipe for producing discretizations sharing the Lax matrix with the 
continuous-time system, so that the discrete-time system belongs to the 
same integrable hierarchy as the underlying continuous-time one [||, f^], 

0. 

In the present paper we want to describe a new application of this scheme 
to the Bogoyavlensky lattices [|lTl, which were given an r-matrix interpre- 



tation in [|I2[. Some of equations derived here appeared previously in the 



literature [0, as certain reductions of the discrete KP equation in the bilin- 



ear form. Our approach enables to get these equations systematically, and, 
moreover, provides automatically the Hamiltonian formulation along with 
the interpolating Hamiltonian flow, as well as the solution in terms of matrix 
factorizations. 



2 Continuous— time Bogoyavlensky lattices 



The Bogoyavlensky lattices were introduced in |[TT|] as three families of inte- 
grable lattice systems depending on integer parameter m > 1 [m > 1 for the 
third one): 

(m m \ 

ak+j - ttk-j , (2.1) 
i=i i=i / 

(m m \ 

n ^k+j - n ^k-j , (2.2) 
3=1 i=l / 

(m m \ m m 

n ^kij - n = n ^kij - n ^k-j- (2.3) 
j=0 j=0 J j=l j=l 

We shall call these systems lattice 1, lattice 2, and lattice 3, respectively. 
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The lattices 1 and 2 serve as generalizations of the famous Volterra lattice, 

ak = ak{ak+i - ak-i), (2.4) 



which is m = 1 special case of both the systems ( |2.1| ), ( |2.2| ). Some special 
case of the lattice 1 was found also independently by Itoh 

The lattice 3 after the change of variables t— and t ^ —t turns 
into 

(m m \ 

n "fc+i ~ n ^k-j , (2.5) 
i=i i=i / 

which serves as a generalization of the so-called modified Volterra lattice, 
the m = 1 particular case of ( p.5|) : 

«fc = a^(afc+i - ak-i). (2.6) 

All these systems may be considered on an infinite lattice (all the sub- 
scripts belong to Z), and admit also periodic finite-dimensional reductions 
(all the subscripts belong to Z/NZ, where N is the number of particles). 
The lattices 1 and 2 admit also finite-dimensional versions with boundary 
conditions of the open-end type: 



for system (|23|) : = for A; < 0, A; > - m + 1; 

for system : = for k < 0, k > N. 

Bogoyavlensky has found also the Lax representations for these systems of 
the form 

f=[T,B], (2.7) 

where for the system ( p.l|) 



T(a, A) = A "^^QkEk^k+m + ^^Ek+i^k, (2.8) 

B{a, A) = ^(o-fc + ttk-i + . . . + ak-m)Ek,k + A*"^^ ^ Ek+m+i,k, (2.9) 
for the system ( p.2|) 

T{a,\) = \ ^^akEk^k+i + ^"'^Ek+m,k, (2.10) 

B{a,X) = —\ " ^ QfcQ.fc+1 ■ ■ ■ o-k+mEk^k+m+i, (2-11) 
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and for the system ( p. 3D 

(2.12) 

B{a, X) = (^k'^o.k+i ■ ■ ■ ak+m-i^k+m^k, (2.13) 

Here for the infinite lattices all the subscripts belong to Z, for the pe- 
riodic cas eall the subscripts belong to Z/iVZ, and for the open-end case 
all the subscripts belong to 1, . . . , A^. Moreover, in the infinite-dimensional 
and open-end cases the dependence on the spectral parameter A becomes 
inessential and may be suppressed by setting A = 1. Below we consider only 
finite lattices. 

All the Bogoyavlensky lattices are Hamiltonian systems. More precisely, 
each system (|2.1| ), (|2.2| ), ( p.3| ) is Hamiltonian with respect to a certain 
quadratic Poisson bracket 

{ofc, ttj} = TTkjakaj, (2.14) 
with a skew-symmetric matrix (vr^j). The corresponding Hamiltonians are: 
H{a) = tr(T'^+^)/(m + 1) = ^ for the systems (g^), 

H{a) = tr(T™'^^)/(m + 1) = X! '^kO.k+i ■ ■ ■ flfc+m-i for the systems 
H{a) = -tr(T--)/m = ^ (^k '<i ■ ■ ■ <^klm for the system ^ 

The Poisson brackets ( p.l4| ), i.e. the matrices (vTfcj), in the context of 
infinite systems were found for the lattice 1 in the original papers by Bo- 
goyavlensky [jll|, and for the lattices 2 and 3 - in Ref. |jl5[. For the finite 



lattices, where some subtleties come out, this was done systematically in [12 



3 Discrete time Bogoyavlensky lattices 

We present now equations of motion of some difference equations which can 
be considered as analogs and approximations to the Bogoyavlensky lattices 
for the case of the discrete time. The "Proposition fc" (A; = 1, 2, 3) deals with 
the " discrete time Bogoyavlensky lattice /c" . We use tilde to denote the time 
shift, so that, for example, Vk = Vk{t + h), if Vk = Vkit). 
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Proposition 1. The system of difference equations 

m m 

^fc 11(1 + hvk-j) =VkW{l + hvk+j) 

j=l 3=1 

admits a Lax representation 

f = L-^TL 

with the matrices 

T(v, A) = A ^ akEk,k+m + ^ X! Ek+i,k, 
L{v, A) = ^ PkEk,k + hX^+' Yl Ek+m+i,k, 

where 

m m 

j=l i=0 

Proposition 2. The system of difference equations 

(m \ I m 

1 + /i n ^k-jj =vk 1 1 + ^ n ^k+j 

admits a Lax representation 

f = UTU-^ 

with the matrices 

T{v, A) = A ^ ^ flfc-Efe^fc+i + A" ^ Ek+rn,k, 
Uiv,X)^I + h\-^-' lkEk,k+m+l, 

where 

(m \ m 

1 + /i n , = n Vk+j. 
3=1 J j=0 

Proposition 3. The system of difference equations 

1 rii \ I in 



admits a Lax representation 

f = L-^TL 

with the Lax matrices 

Tiv, A) = J2 (^kEk,k+i + A-'"-i J2 Ek,k+m+i, (3.10) 

L{v, X) = I + hX^Yl (^kEk+m,k, (3.11) 

where 

1 + ^ n ' = n ^kij- (3-12) 

j=0 / j=0 

Remark 1. Upon change of variables Vk ^ v^^ and h i-^ —h the system 
^ turns into 



Vky-- hY[vk-jj = Vk y- - hY[vk+jj , (3.13) 
which may be considered as a discrete time analog and approximation to 

rsi). 



Remark 2. The equation (p.l|) was found in ||T3[ as a certain reduction 
of the discrete KP equation in the bilinear form. Other equations ( |3.5| ), 
( |3.9|) seem to be new. The equations ( p.l|) and (|3.5|) for m = 1 coincide, as 



they should (Volterra lattice). The Lax representation for this case with the 



matrices (|3.6| ), (|3.7|) was also given in |13], but without any hint on how it 
was obtained. 

In the above formulation these Propositionss may be easily checked by a 
direct computation, but their origin remains hidden. In the following sections 
we shall give a way to derive them systematically, which, as a by-product, will 
unvail an underlying invariant Poisson structure of these discrete systems, 
as well as a role of the auxiliary matrices L, U. This, in turn, will enable 
us to solve the initial value problems for our systems in terms of matrix 
factorizations and to find interpolating Hamiltonian flows. Our construction 
is just a particular case of a general one, applicable, in principle, to every 
system admitting an r-matrix interpretation. The key observation is that the 
Lax matrices ( p.2| ), (|3.6| ), ( p.lO| ) of the discrete time systems formally coincide 
with the corresponding Lax matrices ( |2.8|) , (|2.10|) , ( p. 12]) of the continuous 
time ones. 
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4 Algebraic structure 

of Bogoyavlensky lattices 

In we gave an r-matrix interpretation of the Bogoyavlensky lattices as 
simplest representatives of integrable hierarchies on associative algebras. The 



main results of |T^ may be summarized as follows. 

1) For the open-end case (applies only to the lattices 1 and 2) we set 
g = gl{N). To this algebra there corresponds a group G = GL{N). As a 
linear space, g may be represented as a direct sum of two subspaces, which 
serve also as subalgebras: g = g+ © g_. Here g+ (g_) is a space of all 
lower triangular (resp. strictly upper triangular) by matrices. The 
corresponding subgroups: G+ (G_) is a group of all nondegenerate lower 
triangular A^ by A^ matrices (resp. upper triangular A^ by A^ matrices with 
unities on the diagonal). 

2) For the periodic case (of all lattices 1, 2, 3) g is a certain twisted 
loop algebra over gl{N), namely the algebra of formal semi-infinite Lau- 
rent series T(A) over gl{N), satisfying flT{\)fl~^ = T{ujX), where fl = 
diag(l, tt;, . . . , a;^~^), u = exp{27ii / N) . The corresponding group is the 
twisted loop group G consisting of GL(A^)-valued functions T(A) of the com- 
plex parameter A, regular in CP^\{0, oo} and satisfying flT{\)Q^^ = T{u\). 
Again, as a linear space g = g+ © g_, where for the lattices 1 and 2 g_|_ 
(g_) is a subspace and subalgebra consisting of T(A) containing only non- 
negative (resp. only negative) powers of A, and the case of the lattice 3 
differs in that to which subalgebra do diagonal matrices belong: g_|_ contains 
only positive, and g_ only non-positive powers of A. For the lattices 1 and 
2 the corresponding subgroups G+ and G_ consist of T(A) regular in the 
neighbourhood of A = (resp. regular in the neighbourhood of A = cxd and 
taking the value J in A = oo). For the lattice 3 G+ is formed by T(A) regular 
in the neighbourhood of A = with T(0) = /, and G_ is formed by T(A) 
regular in the neighbourhood of A = oo. 

For both the open-end and periodic cases every T G g admits a unique 
decomposition T = l(T) + u(T), where l(T) G g+, u(T) G g-. Analogously, 
for the both cases every T G G from some neighbourhood of the group 
unity admits a unique factorization T = C(T)U{T), where JC(T) G G+, 
W(T) G G^. 

There hold the following statements. 
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a) For each system ( p.l| ), ( |2.2| ), ( p.3| ) t/iere exists a quadratic r-matrix 
Poisson bracket on g whose Dirac reduction to the corresponding set of ma- 
trices V = {T(a, A)} from ( ^.101 ), or ( |2.12| ), respectively, is given by 

b) Let : g I— > C 6e an invariant function, so that [] d(f{T) = T'S/ip{T) = 
Vip{T)T is covariant under conjugation. Then the Hamiltonian flow on g 
with the Hamiltonian function Lp{T^)/p (here and below p = m + 1 for the 
lattices 1,2, and p = m for the lattice 3) is tangent to V and has the Lax 
form 

t = [TJid^irn)] = - [T.uid^iTn)] . (4.1) 
This flow admits the following solution in terms of the factorization problem 

^tdip{TP{0)) ^ c{t)U{t), C{t) G G+, U{t) G G_ 

(this problem has solutions at least for sufficiently small t): 

Tit) = C-\t)T{0)C{t) =U{t)T{0)U-\t). 

c) Let f : g G be a conjugation covariant function on g. Then the 
difference equation 

f = C-^(^f{TP))T C^fiTP)) = u(^f{TP))TU-^[f{TP)) (4.2) 

defines a Poisson map g g which leaves V invariant, the restriction of 
this map on V being Poisson with respect to the reduced bracket ( p.l4| ) . This 
difference equation admits following solution in terms of the factorization 
problem 

/"(TP(0)) = C{nh)U{nh), C{nh) G G+, U{nh) G G_ 



^Recall that the gradient 'V(p{T) e g of the function 1^9 : g ^ M is defined in the 
open-end case by the relation 



ae 



VXeg; 

e=0 



in the periodic case "tr" should be replaced by "tro", the free term in the Laurent series 
for the trace. 
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(this problem has solutions for a given n at least if /(T(0)) is sufficiently 
close to the group unity I): 

T{nh) = C'\nh)T{0)C{nh) =U{nh)T{<d)U-\nh). 

d) The solutions of the difference equation (|4.2| ) are interpolated by the 
flow ( [4 .11) with the Hamiltonian function (p(TP)/p, where f(T) is defined by 

d^{T) = h~Hog{f{T)). (4.3) 



The statements a),b) explain the Lax equation ( f^.T] ) with the matrices 
(|3)-(|l3D, as for the system we have B{a, A) = /(T'"+i(a, A)), for the 
system ( |2.2| ) we have B{a, A) = —u{T"^^^{a, A)), and for the system ( |2.3| ) we 
have B{a,X) = /(T-"^(a, A)). 



5 A discretization of the 
Bogoyavlensky lattice 1 

We get a correct perspective for the interpretation of the system ( ^TTD (as 
well as the systems ( |3.5| ), ( |3.9| )) if we take an "inverse" view-point. We 
consider the first equation in ( p.4| ) as an implicit definition of the functions 
Vk = Vk{a), rather then the expressions of through Vj. In the open-end 
case the sequence of v^s can be computed even explicitly, term by term, 
starting with Vk = cifc/(l + ^Z]j=i for 1 < A; < m + 1. In particular, for 
m = 1 one has Vk = ak/{l + hvk^i), which implies a nice representation in 
form of a finite continued fraction: 



1 + 



1 + 



ha2 
1 + hai 



In the periodic case the existence of the functions Vk = Vk{a), at least 
for h small enough, follows from the implicit functions theorem. Again, for 



8 



m = 1 we get an expression in the form of an infinite A^-periodic continued 
fraction of the type (^.11) . 

The second equation in ( ^.41 ) may be rewritten as a recurrent relation 
for = /5fc(a). In fact, we have Pk — ha^ = YiJLii^ + hvk-j), so that 
Ofc/ {(3k - hat) = Vk, and finally 

Pk - ha, = n f 1 + , ] ■ (5.2) 

j=i \ Pk~j - tiak-j J 



Conversely, the last formula implies ( |3.4| ), if one sets Vk = ak/{/3k — ha,). 

The formula (|5^ ) may also serve for a successive computation of PkS 
in the open-end case, and in the periodic case it uniquely defines a set of 
f3k — hak, ^ ^ k < N, via the implicit functions theorem. In both cases it is 
easy to see that 

m 

Pk = l + hJ2ak^j + 0{h^). (5.3) 

j=0 

Theorem 1. The quantities (3k defined by ( |5.2| ) serve as coefficients of 
the matrix 

L = C{I + hT^+^) = PkEk,k + h\^+^ Ek+m+i,k. (5.4) 

k k 

The discrete time Lax equation 

f = L-^TL = C-\I + /iT™+^) T £(/ + /iT"^+^) (5.5) 



with the Lax matrix ( |2.8|) generates the following map on M.^ {a}, equivalent 
to (O); 

~ (3k+m fr r\ 

Qk = —z Ofc. (5.6) 

This map is Poisson with respect to the Poisson bracket (|2.14|) corresponding 
to the lattice 1, and is interpolated by the flow with the Hamiltonian function 

^^tr$(T™+i), where ^{i) = h-^ (\og{l + h7])'^ . (5.7) 
m + 1 JO ri 

Proof. The last two statements follow from the results formulated in 
the previous section, provided the first two statements are proved. Suppose 
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for a moment that the £-factor of / + hT^^^ has the form ( |5.4| ). Then the 
evolution equation (|5.5| ), i.e. LT = TL, is equivalent to: 



PkO-k — 0,kPk+m, httk + Pk+m+1 — ha^+m+l + Pk+m- (5. 



This in turn is equivalent to a combination of an evolution equation (|5.6| ) 
with the condition of compatibility of two equations in (|5.8|): 



Pk — hat — J^'^ (Pk+m+i — hak+m+i)- (5.9) 
The last equation is equivalent to the fact that 

m 

Y[i(3k+j - hak+j) 

'^^^^i = const, (5.10) 

n (^k+j 

3=0 

i.e. does not depend on k. We shall prove that the actual value of this 
constant is equal to 1, which is just equivalent to (|5.2| ). 

The inspection of the structure of the matrix T™'^^ for T from ( |2.8| ) 
convinces that the £-factor of J + /iT™"*"^ has in fact the form ( ^.4|) , while 
the W-factor has the form 

m 

i=i k 

The quantities (3k, Ik^ are completely defined by the set of recurrent relations 
following from the definitions: 

m 

Pk + h'jl}l^_, = l + hJ2ak-,, (5.11) 

j=0 

/3fc7i'V/i7En-i = coef. by X'^^"'^'^^ Ek,k+jin.+i) in T^+S l<j<m-l; 

m 

pkiir^ = n (^k+jm- 

j=0 

Now we are in a position to prove that the constant in (|5.10| ) is equal to 

1. 
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Indeed, in the open-end case it is enough to compute from ( p.llD the first 
m + 1 values of namely Pk = ^ + h Y.j=i aj, 1 < k < m + 1, which implies 

nr=V(/3,-K)/nf=i/5, = i. 

In the periodic case we have found only a combinatoric proof based on 
tedious computations. For the sake of simplicity and in order to avoid compli- 
cated notations we present the corresponding argument only in the simplest 
cases m = 1,2. 

In the case m = 1 the defining recurrent relations take the form: 
+ = l + hak + hak-i, Pklk^ = atak+i. 

Excluding 7^^'' from these relations, we get: 

1 _ /3 h /, 

-L — Hk+2 — iT'O'k+2 — iT-O-k+l n • 

Pk 

Replacing the fraction on the right-hand side through its expression following 
from ( p. 91) for m = 1, we get: 

^ _ {.Pk+2 — h ak+2){Pk+i — hak+i) 



which proves the theorem in the case m = 1. 

For m = 2 the defining recurrent relations take the form 

Pk + h'^lklz = 1 + /ittfc + hak-i + hak-2, 

Pklk + ^7fc-3 = O'kak+2 + CtkO-k+l + Ofc-iOfc+i, — 0'kO'k+20'k+4- 

Excluding from these relations we get: 

/3k — hak , , 2 Pk-3 ~ httk- 



Pk+3 — hak+z — h{ak+2 + o-k+i) 3 h h ak+iak-i 



Pk PkPk-3 

According to ( ^.91) for m = 2, this is equivalent to 

1/3 7, 7,/ , -sPk+s - hak+3 2 Pk+3-hak+3 

J- — Pk+3 — ftOfc+S — + a-k+l) n r 



Pk+2 Pk+2Pi 



k-1 



{Pk+3 - hak+3) {Pk+2 - hak+2) , {Pk+3 - hak+3){Pk-i - hak-i) 

nak+v 



Pk+2 Pk+2P, 



k-1 
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Using in the last term once more ( p.9| ) for m = 2, we obtain 
^ _ {Pk+3 — hak+z){(3k+2 — hak+2){Pk+i — huk+i) 

which proves the theorem for m = 2. The pattern of the proof for a general 
m may be seen from these two particular cases. 

6 A discretization of the 
Bogoyavlensky lattice 2 

For the lattice 2 we again consider the first equation in ( p. 81 ) as a defini- 
tion of the functions Vk = Vk{a). In the open-end case we can compute 
these functions succesively, starting with Vk = ak{l + hYl'jZiI\i=i<^i)/i^ + 
^Sj=i n/=i a/) for 1 < A; < m + 1. In the periodic case the implicit functions 
theorem has to be invoked. In particular, for the case m = 1 we obtain the 
same continued fractions expressions as in the previous section. 

The second equation in ( |3.8D may be represented as a recurrent relation 
for 7fc = 7fc(a). Indeed, we have — h'jk-m = Vk, so that 

at - h'jk-m = m • (6-1) 

1 + /i n ("-k-j - hjk-m-j) 

Conversely, the last formula implies p.8|), if one sets Vk = — h'^k-m- 

In the open-end case the formula (|6.1|) serves as a basis for successive 
computation of 7fc's, and in the periodic case it uniquely defines, by the 
implicit function theorem, the quantities ak+m — h'jk, 1 < k < N. In both 
cases there holds the following asymptotic relation: 

m 

-fk = Yl^k+j{l + 0{h)). (6.2) 

i=o 

Theorem 2. The quantities 7^ defined by ( p.lD serve as coefficients of 
the matrix 

U = U{I + hT^+^) = 1 + /iA-(™+i) lkEk,k+m+i- (6.3) 

k 
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The discrete time Lax equation 



T = UTU-^ = U{I + /ir""+^) T U-\I + hT'^+^) (6.4) 



with the Lax matrix ( f^.lOD generates the following map on M^{a}, equivalent 
to 

Ofc — 7 a/c+m+i. lO-C'J 

This map is Poisson with respect to the Poisson bracket ( |2.14| ) corresponding 
to the lattice 2, and is interpolated by the flow with the Hamiltonian function 



Proof. Again, it suffices to prove the ffist two statements. Assuming for 
a moment that the W-factor of the the matrix I + hT"^~^^ for T from ( p.lOp 



has the form (|673|), we see that the evolution equation (|6.4[), i.e. TU = UT, 
is equivalent to 

ak'fk+i = IkCik+m+i, ak + h-fk-m = ak + h'jk. (6.6) 



This in turn is equivalent to a combination of an evolution equation (|6.5| ) 
with the condition of compatibility of two equations in 



'yk 

Ofc - /i7fc_m = {ak+m+i - h-ik+i)- (6.7) 

Ik+i 

The last equation is equivalent to the fact that 

— W{ak+m-j - hjk-j) = const, (6.8) 

i.e. does not depend on k. We shall prove that the actual value of this 
constant is equal to 1, which is equivalent to (|6.1|). 

This time the inspection convinces that the Z//-factor of the matrix / + 
j^j^m+i rp (|2.10| ) must indeed have the form ( p.3|) , while the £-factor 
must have the form 



L = C{I + hT^+^) = (3f^Ek,k + hY, A^-("+i) 

k j=l k 
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where jSj^^ = 1, and other quantities 7^, jS^^^ are completely defined by the 
recurrent relations following from the definitions: 

m 

(^k\k = n«'^+.' (6-9) 

k m— 1 

/3f + /^'4-Ul7.-r.-l = l + h Y: E^l+J, (6.10) 

l=k—m j=0 

f3\^hhl3\^:^Zilk-m-i = coef. by X^'^^^'^ Ek+^(m+i),k in T"^+\ 1 < j < m-1. 



To prove that the constant in (|6.8| ) is equal to 1, in the open-end case 
is enough to compute from (|6.9|), ( |6.1(]| ) the first m + 1 values of 7^, namely 
Ik = nf=oafc+i/(l + ^Etinf="o^a«+j), I < k < m + 1, which implies 

Uf=iiaj+m - /i7i)/7m+i = 1- 

In the periodic case we shall again give the proof only for m = 1,2, leaving 
the tedious calculations for the general case to the reader. For m = 1 the 
defining recurrences (|6.9| ), ( |6.1(]| ) take the form: 

f^k^lk = ctkCik+i, + h'^lk-2 = l + hak + httk-i. 

Excluding from these relations l3f\ we get: 

1 = — (ctfc+i - ^7fe) - ^(ofc-i - /i7fc-2)- 

Ik 

Replacing the last term on the right-hand side through its expression follow- 
ing from (|6.7|) for m = 1, we get: 

. _ (ofc+i - h-fk){ak - h-fk-i) 
Ik 

which proves the theorem for m = 1. 

In the case m = 2 the recurrent relations ( |6.9| ), ( |6.1(JD take the form 

f^k^lk = afcafc+iafc+2, (3^^^ + h'^ pj^lr^-fk-s = l + h{ak-2ak~i + ak-iak + akak+i), 

f^k^ + h^k-3 = ak-i + ak+i + ak+3- 
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Excluding jSjf^ from these relations, we get: 

1 = "^"^"^^ (afc+2 - h-fk) - Kak-2 + afc)(afe_i - h'-fk-z) + h'^'jk-siak-A - h^jk-a) 
Ik 

Using on the right-hand side repeatedly ( |6.7|) for m = 2, we can rewrite it 



as 



1 = [ak+2-hlkj [ak+2-hlk>^ [ak+2-hlk> 

Ik Ik Ik 

= — (afc+2 - h'jk){ak+i - h'jk-i) - - /i7fc)(afe-2 - h'^k-i)- 

Ik Ik 

Using in the last term once more ( |6.7D for m = 2, we get 

1 = —{ak+2 - h-fk){ak+i - h-fk-i){ak - /i7fc-2), 
which finishes the proof for m = 2. 

7 A discretization of the 
Bogoyavlensky lattice 3 

For the lattice 3 we again define the functions Vk = Vk{a) by means of the 
first equation in (|3.12 ), which is justified by the implicit function theorem (as 



opposed to the lattices 1,2, this time an open-end reduction is not admissible, 
so that only the periodic case needs to be considered). In particular, for 
m = 1 we have Vk = ak — h/vk-i, which leads to the expression in terms of 
an infinite A^-periodic continued fraction: 

h 



Vk = ak 



dk-r 



h 



h 



Vk 



The second equation in ( p.8| ) implies Uk — hak-m = Vk, and hence 

I 

=0 ^k+j — hak+j 



m—l -1 

«^ = n . • (7-1) 
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Conversely, the last formula implies ( |3.12| ), if one defines Vk = ak — hak-m- 

The formula ( [7.1|) defines, by the implicit function theorem, the set of 
quantities ak, 1 < k < N, satisfying 

m— 1 

«fc= n%i(i+ow). (7.2) 

j=0 

Theorem 3. The quantities ak defined by ( |7.1D serve as coefficients of 
the matrix 

L = C{I + hT-"') =I + hX^Y. <^kEk+m,k. (7.3) 

k 

The discrete time Lax equation 

f = L~^TL = C-\I + hT-"^) T C{I + hT-"^) (7.4) 

with the Lax matrix (p.l2| ) generates the following map on R^{a}, equivalent 
to (PI).- 

~ (^k hak—m (rj r\ 

Ofc — 7 Ctfc+m.. V -^j 

ak+m — nak 

This map is Poisson with respect to the Poisson bracket (|2.14|) corresponding 
to the lattice 3, and is interpolated by the flow with the Hamiltonian function 

-— tr<l>(T-™), where $(0 = h'^ log(l + /ir/)^. 
m Jo rj 

Proof. Again, it suffices to prove the first two statements. Assuming for 
a moment that the £-factor of the the matrix / + hT~^ for T from ( |^.12| ) 
has the form (|7.3| ), we see that the evolution equation (|7.4|) , i.e. LT = TL, 
is equivalent to 

akak = ak+mOk+i, ak + hak-m = ak + hak+i- (7.6) 

This in turn is equivalent to a combination of an evolution equation ( |7.5| ) 
with the condition of compatibility of two equations in (|7.6| ): 



ttk - hak-m = ^^^^{ak+m - hak). (7.7) 
ak 



The last equation is equivalent to the fact that 



m— 1 

ak n ('^fc+i " hak+j-m) = const, (7.8) 
i=o 
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i.e. does not depend on k. We shall prove that the actual value of this 
constant is equal to 1, which is equivalent to (|7. 1|) . 

To compute the £-factor of the matrix / + hT~'^ for T from ( ^.12| ) , we 
notice, first, that = CD~^, where 

C = \ ^ a^^Ek+i,k, D = I + A""" a^l^Ek,k+m- 

k k 

Further, notice that the ^-factor of any matrix is not changed under the right 
multiplication by the factor from G_. We multiply the matrix / + hT~^ = 
I + (CD~^)^ from the right by (DC^^)"^C"^. To see that this matrix belongs 
to G_, notice that it is equal to DDi . . . Dm-i, where Dj = C~^DC^ = 
/ + A~™ J2 dk^Ek^k+m £ G_. For the further reference we give here an explicit 
formula 

J— 1 j 

dk = n ^^+1 n ^k+m+i- 

1=0 1=0 

So we get 

C{I + /iT-™) = CiDDi . . . An-i + hC""), 

and an inspection of this formula convinces that this factor must indeed be 
of the form (|7.3|) , while 

m 

U{DD^ . . . D^_, + hC^) = J2 A-^"^ E Pk^Ek,k+jm 

j=0 k 

Here the quantities at, 13): are completely defined by the recurrent relations 
following from the definitions: 

m— 1 

c^kPT = n ^klr (7.9) 

j=0 

+ ha,.^/3j^l^ = 1, (7.10) 
tf'^ + = coef. by X-'^'Ek^k+jm in DD,... 1 < j < m. 

(In the last equation for j = m one must set (3^^^^^ = 0, which leads to 

Again, we shall prove that the constant in ( [7.8| ) is equal to 1, only for the 
two simplest cases m = 1,2, leaving the calculations for the general case to 
the reader. 
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For m = 1 the defining recurrences ( [7.9| ), ( [7 .101 ) read: 

„ /5(0) _ ^-1 /?{0) I /.^ oW _ 1 _ „-l 

OfcPfc — O-k ■> Pk + iT'O^k-lPk-l — J- 5 Pfc — 

Excluding from these relations (3^\ l3j^\ we get: 

— = 1 - /ittfc or — = ak- hak-i, 

Oik ak 

which proves the theorem for m = 1. 

In the case m = 2 the recurrent relations ( [7 .91) , (|7.1CI| ) take the form 

o^kP^k^ = ajfc^ctfc+i, Pj^^ + hak-2l3kl2 = 

f^i '* + hak-2(^k-2 — + (^kO'k+2^k+3^ f^k ^ — ^k+4^k+5- 

Excluding jS^^^ from these relations, we get: 



'^k ^Q-fc+i 

Oik 

or 



1 - hak-2{ak ^ + afc-2a.fc ^ctfc+i - hak-^a^ ^ctfc+i)) 



— = ak+i{ak - /iafc_2) - hak-2{ak-2 - hak-4). 

Oik 

Using in the last term on the right-hand side ( [7.7| ) for m = 2, we can rewrite 
the last expression as 

— = afc+i(afc - hak-^2) - hak-i{ak - hak-2) = {ak+i - hak-i){ak - hak-^2)- 

Oik 

This finishes the proof for m = 2. Again, we hope that the pattern of 
the general proof is clear from these two simple cases. It would be highly 
desirable to find a less computational proof for the periodic case of all three 
lattices. 



8 Conclusion 

A new application of a general scheme for producing integrable discretizations 
for integrable Hamiltonian fiows is described in the present paper. Advan- 
tages of this approach are rather obviuos: it is, in principle, applicable in 
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a standartized way to every system admitting an r-matrix formulation, at 
least with a constant r-matrix satisfying the modified Yang-Baxter equa- 
tion. We shall demonstrate elsewhere that the discrete time systems from 
0) with dynamical r-matrices may be also included into this framework. 
We hope also to report on numerous further applications of this approach in 
the future. 

The drawback of this scheme is also obvious to any expert in this field. 
Namely, some of the most beautiful discretizations do not live on the same 
r-matrix orbits as their continuous time counterparts 0, @|, p, and 
there seems to exist no way of a priori identifying the correct orbit for nice 
discretizations. However, we hope that continuing to collect examples will 
someday bring some light to this intriguing problem. 

The research of the author is financially supported by the DFG (Deutsche 
Forschungsgemeinshaft) . 
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